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Abstract 

Searching for the integrable structures of supersymmetric gauge theories and topolog- 
ical strings, we study melting crystal, which is known as random plane partition, from the 
viewpoint of integrable systems. We show that a series of partition functions of melting 
crystals gives rise to a tau function of the one-dimensional Toda hierarchy, where the mod- 
els are defined by adding suitable potentials, endowed with a series of coupling constants, 
to the standard statistical weight. These potentials can be converted to a commutative 
sub-algebra of quantum torus Lie algebra. This perspective reveals a remarkable connec- 
tion between random plane partition and quantum torus Lie algebra, and substantially 
enables to prove the statement. Based on the result, we briefly argue the integrable struc- 
tures of five-dimensional N = 1 supersymmetric gauge theories and ^-model topological 
strings. The aforementioned potentials correspond to gauge theory observables analogous 
to the Wilson loops, and thereby the partition functions are translated in the gauge the- 
ory to generating functions of their correlators. In topological strings, we particularly 
comment on a possibility of topology change caused by condensation of these observables, 
giving a simple example. 
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1 Introduction 



An unanticipated but very exciting connection between the statistical mechanical problem of 
melting crystal, known as random plane partition, and A-model topological strings has been 
revealed [1] , based on the topological vertex [21 [3] . The topological vertex is a diagrammatical 
method which enables to compute all genus topological A-model string amplitudes for a certain 
class of local geometries. 

We can image the positive octant Z>q C M'^ occupied by unit cubes as the neighborhood of 
a corner of the crystal by putting unit cubes on the lattice points in the octant. The frozen 
crystal occupies the positive octant Z>q. As the crystal melts, we remove atoms from the corner. 
We identify the configuration of crystal melting as the configuration of plane partition or the 
three-dimensional Young diagrams, as depicted in Figure [H Removing each atom contributes 
the factor q = e~T to the Boltzmann weight of the configuration, where T is the temperature 
and /i is the chemical potential. Heating up the crystal leads to melting of it. 




Figure 1: The corner of the melting crystal and the corresponding plane partition (the three- 
dimensional Young diagram). 



Random plane partition also has a significant relation with five-dimensional Af = 1 super- 
symmetric gauge theories. Nekrasov's formula [U [5] for five-dimensional A/" = 1 supersymmetric 
SU{N) Yang-Mills theory can be retrieved from the partition function of random plane parti- 
tion [6] , where the model is interpreted as a g-deformed random partition. All genus topological 
A-model string amplitude for the local SU{N) geometry is evaluated by the topological vertex 
and reproduces Nekrasov's formula for A/" = 2 SU{N) gauge theory [7, 8], as predicted in the 
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geometric engineering. 

It is shown in [5j that the Seiberg-Witten solutions [9J of four-dimensional M = 2 super- 
symmetric gauge theories emerge through random partition, where Nekrasov's functions for 
four-dimensional M = 2 supersymmetric gauge theories are understood as the partition func- 
tions of random partition. The integrable structure of random partition is elucidated in [10], and 
thereby the integrability of correlation functions among certain observables in four-dimensional 
Af = 2 supersymmetric gauge theories is explained. 

Motivated by the these results, we study in this article the integrable structure of random 
plane partition in order to search for the integrable structures of supersymmetric gauge theories 
and topological strings. 

A partition A = (Ai, A2, ■ ■ ■ ) is a sequence of non- negative integers satisfying Aj > Aj+i for all 
i > 1. Partitions are identified with the Young diagrams. The size is defined by |A| = J2i>i -^j; 
which is the total number of boxes of the diagram. A plane partition vr is an array of non- 
negative integers 



satisfying vTjj > iVi^ij and TCij > vTjj+i for all i,j > 1. Plane partitions are identified with the 
three-dimensional Young diagrams. The three-dimensional diagram vr is a set of unit cubes 
such that TTjj cubes are stacked vertically on each (z, j)-element of vr. The size is defined by 
|7r| = Ylij>i '^iji which is the total number of cubes of the diagram. Diagonal slices of a plane 
partition vr become partitions, as depicted in Figure 2. Let 7r(m) denote the partition along 
the m-th diagonal slice, where m G Z. In particular, 7r(0) = (vth, 7122, ■ ■ ■ ) is the main diagonal 
partition. This series of partitions satisfies the condition 



TTii 7ri2 TTis 



7r2l 7r22 vr23 



(1.1) 



TTsi Vr32 7133 



< 7r(-2) -< 7r(-l) ^ 7r(0) ^ 7r(l) ^ 7r(2) > , 



(1.2) 



where jj y u means the interlace relation between two partitions /x, u. 




(1.3) 
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(a) 



(b) 



Figure 2: Plane partition (The three-dimensional Young diagram) (a) and the corresponding 
sequence of partitions (the two-dimensional Young diagrams) (b) . 



A statistical model of plane partitions is introduced by the following partition function. 

Z ^ ^gl-l, (1.4) 

TT 

where the sum is over all plane partitions. The parameter q is indeterminate satisfying < < 1 
and play a role of chemical potential (the energy of the removal of an atom from the crystal). 
Summing over plane partitions, we can obtain 

oo ^ 

E = n 71^- (1-5) 

TT n=l ^ ' 

The partition function of the model is the generating function of plane partitions, known as 
the McMahon function. 



1.1 The models 

We introduce a series of potentials for partitions and put the main diagonal partition of tt 
in these potentials. For the later convenience, we introduce them as functions on charged 
partitions, which are partitions paired with integers. Let $fe (/c = 1, 2, ■ ■ ■ ) be the following 
functions on charged partitions.® 

oo oo 
i=l i=l 
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where (A,p) denotes a charged partition. Actually, the right hand side of this formula becomes 
a finite sum by cancellation of terms between the two sums. More precisely, 

M>^,P) = - g^(P-+i)) + q'—^ . (1.7) 

i=l 

With each fixed value of p, these provide a series of potentials for partitions. These functions 
have been exploited in [TU] from the four-dimensional gauge theory viewpoint, with q or g'^ being 
replaced by a generating spectral parameter. Introducing the coupling constants t = {ti,t2, ■ ■ ■), 
we write their combination as 

^it;p){X)=J2'^M^,p). (1.8) 
k=l 

The partition function of random plane partition whose main diagonal partition is in the po- 
tential ( II. Sp is defined by 

Zp(t) = ^gl"le*(*^-)W°». (1.9) 

The model has an interpretation as a g-deformed random partition. To see this, note 
that, by virtue of the interlacing relations (11.21) . the two series of partitions {7r(m)}^_p and 
|7r(— m)}^_jj represent a pair T, T' of semi-standard Young tableaux of shape vr(0), in which 
the part of the m-th skew Young diagrams 7r(±m)/7r(±(m + 1)) is filled with m + 1. The 
partition function can be thereby reorganized to a sum over the Young diagram A = 7r(0) and 
the pair T, T' of semi-standard Young tableaux of shape A as 

Z.it)=Y: E lVe'^-^^^'\ (1.10) 

A r,T' : shape A 
oo oo 

where = H g^^+^^^l'^^^l-l'^^+^^D and q^' = JJ q{m+l)(M-m)\-\ni^m-i)\) _ rj.^^ pg^^^j^^j 

m=0 m=0 

over the semi-standard tableaux gives the Schur function sx{qP) = s\{xi,X2, ■ ■ ■) specialized to 

Xi = q 2 : 

E E 9^' = ^>^in (1-11) 

T: shape A T': shape A 

Therefore the partition function can be eventually expressed as 

Z,(t) = E ^A(gO'. (1.12) 

A 
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The representation of sx{q^) in terms of the hook polynomial [Ilj allows us to write it further 
as 

where h[s) denotes the hook length of the box s G A. 



1.2 Main result 

In this article, we study these models of random plane partition from the viewpoint of integrable 
systems. Our main result is that the following series of the partition functions is a tau function 
of an integrable hierarchy: 

T{t; p) = g^P(P+i)(2p+i)^^(t) ^ p^z. (1.14) 

More precisely, we show in the text that T{t;p) is a tau function of the one-dimensional Toda 
hierarchy, where the coupling constants t = {ti,t2,---) are interpreted as a single series of 
time variables of the one-dimensional Toda hierarchy. In particular, r(t; p) is shown to have a 
representation by two-dimensional free fermions or free bosons as 

T{t;p) = e^"=' (p| e^^"=i(-)''*'='^'= e^^"=i("^''**'^-* |p) . (1.15) 

In the right hand side of this formula, is an element of GL{oo) and J±k denote the modes of 
the U{1) current associated with the complex fermions. Thus T{t]p) is a matrix element, taken 
between the Dirac sea of U{1) charge p, of an element of GL{oo). The exponential operators 
on both hand sides of gj, are generators of the commuting flows of the one-dimensional Toda 
hierarchy. 

Once it is known that the partition function relates with the tau function, it becomes 
amenable to obtain an infinite set of non-linear differential equations that the partition function 
obeys. This is because tau functions of an integrable hierarchy satisfy an infinite set of non- 
linear differential equations. These non-linear differential equations are all encoded in the 
bilinear identity. It follows from the above that the partition functions satisfy the following 
bilinear identities in the one- dimensional Toda hierarchy. 
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= q{p+p'+i)(p-p'+i) ^ ±.^P-P'e-h^.>lii^-tD^-''Zp+,{t+[z])Zp,_^{t' - [z]) , (1.16) 

where p, p' are arbitrary. The integral on the left hand side means taking the residue at 2; = oo 
and multiply it by —1; the integral on the right hand side is understood to be the residue at 

2 3 _ 

z = 0. We also use notations like t ± = (ti ± ^2 ± ^3 -'- T' ' ' ' )• Towers of non-linear 
differential equations are obtained from fll.l6p as the coefficients of the Taylor expansions along 
the diagonal t = t', that is, the coefficients of the expansions of the bilinear identities in the 
variables tk — t'/^. For instance, let p = p' . The first equation one gets in the tower is nothing 
but the Toda equation in Hirota's bilinear form. 

D2Zp-Zp = g2f+iZp+iVi, (1-17) 

where D denotes Hirota's derivative that is defined by D.j.]{x) ■ g[x) = lim{dx — dy)f{x)g{y). 

The partition functions also give rise to a solution of the modified KP hierarchy. The 
corresponding bilinear identities are read as 

/ ^,p-p' e^.>.i'u-t',)^' Z,{t _ [;,-i])Zp,(t' + [z-^]) = , (1.18) 

where p > p'. These identities include in the towers the modified KP equations as well as the 
KP equation. For instance, let p' = p — 1. The first equation one gets in this tower is 

{Dl-D,,)Z,{t)-Z^^,{t) = 0. (1.19) 

This is the first equation of the modified KP hierarchy. 

1.3 Transfer matrix approach 

The main tool we use in the text is the transfer matrix formulation of random plane partition. 
The hamiltonian picture is hinted from the interlace relations (11. 2p . which state that plane 
partitions are certain evolutions of partitions by the discretized time m. In particular, the 
transfer matrix formulation [12j makes it possible to express the partition function (11.91) in 
terms of two-dimensional conformal field theory {2d free fermions). 

Let ipi^z) = J2m&'^rnZ~"^~^ and ip*{z) = XlmGZ^m^"'" ^c complex fermions with the anti- 
commutation relations, {ipm, i^n} = ^m+n,o and {ipm, i^n} = {i^m, i^n} = 0- The Noether current 
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of the f/(l) rotation is given by 

J{z) = : ^{z)r{z) : = z-'^-'Jm , (1.20) 

where : : denotes the normal ordering of fermions that is defined by 

iIj{z)^*{w) =:tp{z)tp*{w) : +^—, \z\ > \w\ . (1.21) 

z — w 

It is well known that partitions are realized as states of the fermion Fock space. In particular, 
charged partitions are realized by the states of the same U{1) charges. For a charged partition 
(A,p), the corresponding state is 

oo 

1=1 

where \p) denotes the Dirac sea having the U{1) charge p, and is defined by the conditions 

i^mlp) = for Vm > —p, 

^^\p) = forVm>p + l. (1.23) 

Using the realization (11.221) . it can be seen that the function (11. 6p corresponds to the fol- 
lowing fermion bilinear operator. 

H, ^ 5^ g'™ : V'-n^C : • (1-24) 
Actually, the potential function (11.61) is reproduced as 

Hk\X;p) = ^kiX,p)\X;p). (1.25) 
These operators are commutative. The following combination reproduces the potential (II. 8p . 

oo 

H{t) = Y,hHk. (1.26) 

k=l 

The transfer matrices [12j are vertex operators of the following forms. 

+ 00 

r+(m) = exp(^-g-'=(-+5)j,), (1.27) 

k=l 
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+ 00 

r-M = exp(^-g^'("+^)j_,), (1.28) 

k=l 

where J±k are the modes of the U{1) current. The matrix elements between partitions of 
different charges always vanishes, while those between partitions of the same charge are 

-{m+i)(|^t|-|A|) A -< U 

(A;p|r+(m)|/x;p) = { ^ ^ (1.29) 

otherwise . 

{fi;p\T4m)\X;p) = { ^ ^ (1.30) 

otherwise . 

By comparing these formulas with the interlace relations (11. 2p . we see that T±{m) describes 
the evolutions of partitions. More precisely, the evolution at a negative time m < — 1 is given 
by r_|_(m), while the evolution at a nonnegative time m > is by r_(m). 

Taking the hamiltonian picture of plane partitions, the partition function (11.91) can be repro- 
duced in the transfer matrix formulation. Actually, following the same steps as we translated 
the partition function to the g-deformed random partition (1.10), but using the transfer matrices 
in place of the Schur functions, the partition function (11. 9p is eventually expressed as 

Z,{t) = (p|G+e^WG_b), (1.31) 

where G± are the propagators those are responsible respectively to the negative time evolutions 
and the nonnegative time evolutions of partitions. These are the operators given by the following 
infinite products. 

-1 

G+ ^ n r+M, (1-32) 

m=— oo 

+ 00 

G_ = JJr_(m). (1.33) 

m=0 

1.4 Quantum torus Lie algebra and random plane partition 

Starting from the expression (ll.3ip of the partition function (11.90 . we prove in the text that 
the series of the partition functions (I1.14p is a tau function of the one-dimensional Toda hier- 
archy. Remarkable relations between random plane partition and quantum torus Lie algebra 
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are revealed in the course of the proof. Throughout the text, we take the perspective that such 
a quantum Lie algebra is a hidden symmetry of random plane partition. 

We realize the quantum torus Lie algebra in terms of the complex fermions. Using this 
realization, we can regard the operators Hk as a commutative sub-algebra of the quantum 
torus Lie algebra. The adjoint actions of the propagators G± on the Lie algebra generate 
automorphisms of the algebra. By taking advantage of such automorphisms, we provide a 
proof of the statement. Actually, among such automorphisms, we pay special attention to the 
shift symmetry that is the automorphism generated by the adjoint action of the product G-G^ 
or equivalently G^G^. By utilizing this symmetry, we can eventually express the partition 
function in the form (ll.lSp . In particular, the element of GL{oo) in the formula (ll.lSp is given 
by 

g. ^ , (1.34) 

where is a generator of the PFoo-algebra of the following form: 

W ^ J^'^'-'^^-mrm-- ■ (1-35) 

Using this formula, we finally confirm the statement by showing that actually realizes a 
solution of the one-dimensional Toda hierarchy. 

Taking the formula fll.lSp . Virasoro/W-constraints on the partition function (11. 9p and the 
tau function (I1.14p can be obtained from the transformation of the PVoo-algebra by the adjoint 
action of g^. In the same way, the transformation of the quantum torus Lie algebra by the 
adjoint action of g^ gives rise to quantum torus analogues of the Virasoro/W-constraints on the 
partition function (II. 9p and the tau function (11.140 . As we argue subsequently, such quantum 
torus analogues of the Virasoro/W-constraints are also obtainable in five- dimensional A/" = 1 
supersymmetric gauge theories and certain topological string amplitudes. These are reported 
in a separate publication [13]. 

1.5 Integrability of 5d J\f = 1 SUSY gauge theories from random 
plane partition 

Random plane partition has a significant relation with five- dimensional A/" = 1 supersymmetric 
gauge theories P, [IH [151 [IS] ■ Our study of the integrable structure of random plane partition 
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is motivated by a quest for the integrable structure of five- dimensional J\f = 1 supersymmetric 
gauge theories and topological strings. 



1.5.1 Random plane partition and 5d J\f = 1 SUSY gauge theories 

Nekrasov's functions for five- dimensional SU (N) gauge theories [H are interpreted as par- 
tition functions of random plane partition [6]. Actually, the original partition function fll.4p 
reproduces Nekrasov's function for five-dimensional U{1) gauge theory, by adding a simple 
chemical potential Q^'^^^^^ to the statistical weight. In the transfer matrix approach, that par- 
tition function becomes 

^5di/(i) _ (0|G+Q^«G_|0), (1.36) 
where Lq = y m : ip^rn'^m '■■ The above matrix element is easily computed and leads 



+ 00 ^ 

11 (1 _ Qqn^n 

Indeterminates q, Q in the right hand side of this formula must be interpreted in terms of 
the gauge theory parameters to reproduce Nekrasov's function for the five-dimensional gauge 
theory. The gauge theory lives on x S"^, where R denotes the radius of the circle, and has the 
dynamical scale A. The indeterminates are identified with these parameters by the relations 

q = e'^^, g = (i?A)2. (1.38) 

By shrinking the circle to a point, from Nekrasov's functions for five-dimensional gauge 
theories, one obtains the four- dimensional versions. Considering the relations (11.38^ . this in- 
dicates that the partition function fll.36p is a g-analogue of the four- dimensional version. To 
see this, we note that the four- dimensional limit of the right hand side of (11.361) is obtained by 
employing the dressed propagators, Q~^^'^G-f-Q^^" and Q^^°G-Q~^^° . Actually, taking the 
relations (11.381) . these dressed operators become nonsingular at the limit R ^ 0, and eventually 
give the following operators. 



lim g-5^«G, g^^o = e^-^' , (1.39) 

lim Q^^^G-Q-^^" = A^o e^-^-i A-^° . (1.40) 



lim 
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By using these formulas, one obtains 

lim(0|G+Q^°G_|0) = lim(0|Q-^^oG+Q^^°Q^^°G_Q-^^»|0) 

R — R — >^ 

= (0| e^^^ A^^" e^^-i |0) . (1.41) 

The right hand side of this formula is nothing but Nekrasov's function for four-dimensional 
U{1) gauge theory. 

1.5.2 Integrable structure of 5d J\f = 1 SUSY gauge theories 

The operator has a counterpart in five-dimensional gauge theories. It corresponds to the 
Wilson loop operator encircling the circle k times |17j . 

Ok = Tr ^Pef^'^^+^'^y , (1.42) 

where and (p denote respectively the fifth component of the gauge field and the real scalar 
field in the vector multiplet. Generating function of the correlators among these observables 
becomes thereby the following analogue of (11.311) . 

^MC/(i)(^) = (p| e^W g^°G_ \p) . (1.43) 

In the same way as we stated on (11.141) . the series of the generating functions (11.431) also gives 
rise to a tau function of the one-dimensional Toda hierarchy. In particular, the tau function 
has the following expression analogous to (ll.lSp . 

.1 ^ g fe=iTr7^ (ple^Er^iH'tfcA gMi/(i) giEr=i(-)"tfe./-fc |p) ^ (1.44) 

where g*'^^''^'' is the element of GL(oo) given by 

gMC/(i) ^ (G'_G+) (G_G+) . (1.45) 

We note that this formula is easily generalized to the SU{N) gauge theory, where the series 
of the corresponding generating functions of the correlators becomes a tau function of the 
one-dimensional Toda hierarchy. 

In four-dimensions, such an integrable structure of A/" = 2 supersymmetric gauge theories 
has been found out [IDl [IE] among the generating functions of the higher Casimir operators 

11 



TT(f)^, where (p is the complex scalar in the vector multiplet, corresponding to A4 + i(p in five- 
dimensional theories. One might expect that the tau function (11.441) is a g-analogue of the 
tau function of the four- dimensional theory, just like the partition function fll.36p is the q- 
analogue. However, this is not the case. Relation between these two integrable structures is 
not straightforward. The subtlety can be found, for instance, in that the tau function fll.44p 
becomes trivial at the four- dimensional limit, owing to the degeneration of all the observables 
(11.421) . Actually, when R is nearly zero, the generating function (11.431) behaves as 

^(^)(t) ~ i?p(p+i)ePJ^^=i*^ (0| e^^' A^^o e^-^"^ |0) . (1.46) 

1.6 Integrability of topological string amplitudes from random plane 
partition 

In addition to the gauge theory interpretation, the partition function (I1.36P has an interpre- 
tation as all genus A-model topological string amplitude on C © 0{—2) CP^ It is a 
non-compact toric Calabi-Yau threefold, often called a local geometry. The toric description is 
given by a fan A consisting of rational cones of dimensions < 3 on M^. The topological vertex 
[21 [3] is a diagrammatical method to compute all genus A-model topological string amplitudes 
for such local geometries. The diagram can be drawn from a polyhedron which is obtained 
by taking duals to the cones in M^. For the local geometry O © 0{—2) CP^, the relevant 
diagram is depicted in Figure [31 The topological vertex computation based on the diagram 
gives all genus A-model topological string amplitude on (9 © 0{—2) — > CP^ as 

+00 ^ 

-^string " Y\. (1 _ (.-a^-ngsAn ' (1-47) 
n=l ^ ' 

where a denotes the Kahler volume of the base CP^, and gst is the string coupling constant. 
Comparing two formulas (11.371) and (ll.47p . one sees that the partition function (I1.36P becomes 
all genus A- model topological string amplitude on O (B 0{—2) CP^, by the following identi- 
fication of the parameters. 

q = e-^'\ Q = e-\ (1.48) 
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Figure 3: The diagram for O © C(-2) CP\ 



As is the case of the five- dimensional gauge theory, we generahze the amphtude including 
the operators in it, anticipating that a wound Euclidean brane along the M-theory circle 
corresponds to the observable Ok of the gauge theories. We physically conjecture such a 
generalization by 

AZfSt'\^-^P) = (P\ G+e^^'^ Q'^'G. \p) , (1.49) 

where the right hand side of this equation is same as the formula fll.43p but with the different 
interpretation fll.48p . The series of the generating functions (ll.49p gives the same tau function 
as is the case of the five-dimensional U{1) gauge theory. 

It seems rather rare that the one-dimensional Toda hierarchy shows up as an integrable 
structure of topological strings. One of such rare cases is the topological sigma model (in other 
words, the Gromov-Witten invariants) of CP^ [HI [20l [211 [221 [231 [2l]. On the other hand, 
the "relative" or "equivariant" versions of those Gromov-Witten invariants have a different 
integrable structure, namely, the two-dimensional Toda hierarchy [231 l25]. It is remarkable 
that substantially the same quantum torus Lie algebra is used in the work of Okounkov and 
Pandharipande [25], but we have been unable to see whether there is a deep connection with 
our work. To obtain the generators Vm^ of our quantum torus Lie algebra, one has to specialize 
the parameter z of Okounkov and Pandharipande's operators Sm{z) to = q^; such powers of 
q apparently play no role in the work of Okounkov and Pandharipande. 

The two-dimensional Toda hierarchy is also known to arise in the generating function r(x, x) 
= J2x 11 sx{x)s^{x)cxfj,, of the two-legged topological vertex cx^, [26]. By changing variables from 
X, X to Tk = ^ J2i ^fc = ~i Si ^i^y ^) coincides with the special part r'^'^°'^{T, T, 0) of a 
^We thank Y. Hyakutake for suggesting such a possibihty to us. 
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tau function r^^°'^"(T, T,p) of the two-dimensional Toda hierarchy. As Zhou pointed out, the 
tau function t'^^°'^"-{T,T,p) has a fermionic representation in terms of an element g of GL{oo) 
of the form 

g = q'^G+G.q^, (1.50) 

where K = ^^(m — ^)^ : ip-m'^^ :. Thus the building blocks of this tau function are similar 
to the foregoing partition functions of U{1) gauge theory. The difference between K and W is 
almost negligible, and can be absorbed by rescaling of T^^Tk and an exponential prefactor. 



1.6.1 Emergence of geometry from condensation 

It is amazing to comment on a possibility of emergence of another geometry from the conden- 
sation in a local geometry. 

Consider the generating function fll.49p for the case of p = 0, and write it simply as 
A^^^g '^\t) = Ji^^^g~'^\t\p = 0). As explained in the beginning of Section 4, ( 11.49^ has 
another representation of the following form. 

AZ^Zl-'\t) = e^^- ^ (01 eS-(-)^'^^^ gf-W 10) . (1.51) 

This becomes a tau function of the KP hierarchy. 

The generating function has a factorization in terms of the skew Schur functions. Ma- 
trix element in the right hand side of formula f ll.5ip is factored by plugging the unity 1 = 
J2p£zJ2x P) {^'i P\ so that it divides g^'^^''^'' into two parts. Owing to the charge conserva- 
tion, the sum over p truncates to p = 0. Writing |A; 0) simply as |A), we can obtain 

(0| e^*°°=i(-)'**-^^- gf ^(1) |0) = J2 (0| e^^-i(-)'*^^* |A) x (A| G+g^«G_ |0) . (1.52) 

A 

Matrix elements in the right hand side of this equation are expressed in terms of the skew Schur 
functions as follows. 



(0|eEr^i(-)'=*^-^'= gfG_|A) = $^g^+^s^(x)s^/A(gO, (1-53) 
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where variables are changed from t to x hj tk = ^ Ylii^'^i)^- Combining these two formulas, the 
right hand side of formula fll.Sip is expressed in terms of the skew Schur functions. Therefore, 
the factorization of the generating function, normalized by the original A-model topological 
string amplitude, can be written as 



.OffiO(-2). X fc 
^string l^J T,^-, ^ 

•^string V"/ 



e 




-.'^ V Q\^\s,iqP) { V 9^^+^ s^{x)s^/,{qP) \ . (1.55) 



^string 

We examine the condensation by choosing the coupling constants t at certain values. In 
particular, we take the following values. 

Owing to the identification fll.48p . t* behave ~ gj^^ when gst is nearly zero and therefore a 
possible condensation becomes nonperturbative in HA superstrings. However, one can compute 
the right hand side of formula (11.551) . The computation yields eventually the amplitude as 



.0®0(-2)/ 
•^string v'' 

"^string \ J n=l 



,0®0{-2) 



JJ(l_gg-+|)". (1.57) 



The right hand side of this formula coincides with all genus A-model topological string ampli- 
tude on the resolved conifold 0{—l) © — » CP^, where the Kahler volume of the base 
is a + ^Qst. 

Emergence of the resolved conifold in (11.571) seems mysterious. However, this can be ex- 
plained as follows. Let us consider the local geometry of coupled conifolds. Coupled conifolds 
can be obtained by patching torically together the foregoing local geometry and C^. The 
diagram of coupled conifolds is depicted in Figure HI where (5i,2 are the Kahler parameters 
attached to the internal edges of the diagram. Each internal edges corresponds to CP^. The 
Kahler parameters are given by Q12 = e"'^^'^, where 2 denote the Kahler volumes of the 
corresponding CP^s. Based on this diagram, the topological vertex computation yields the all 
genus ^-model topological string amplitude as 



A two conifolds 

"^string ~ TT+OO/'-' ^ ' l^i.OOJ 
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Figure 4: The diagram for two conifolds. 



It is remarkable that the topological string amplitude fll.SSp appears, by tuning the Kahler 
volumes ai^2, as a building block of the generating function (11.511) att = t*. Actually, the matrix 
element (0| e^'==i*^~^''*'="^'= g*'^^''^'' |0) is evaluated by using formulas (ll.53l) . (ll.54l) and becomes 
eventually as follows. 

(0|eS^=i(-)'*'=^^ gf^W|0)[_^^ = A'ZZT^°"'\ (1-59) 

where the Kahler volumes ai,2 in the right hand side of this formula are 

1 

ai = a, 02 = -gst- (1.60) 

This formula indicates that the condensation (I1.56P changes the original geometry into two- 
conifolds with two-cycles having the Kahler volumes (11.601) . and that the ratio (11.571) counts 
the worldsheet instantons wrapping the two-cycles simultaneously. The further issue will be 
reported in [13]. 



Organization of the article 

The purpose of this article is to show that the series of the partition functions (I1.14p is a tau 
function of one-dimensional Toda hierarchy. We start Section 2 with giving the realization 
of quantum torus Lie algebra by using the complex fermions. In Section 3, we argue the 
automorphisms of the Lie algebra generated by the adjoint actions of G±. By using such 
automorphisms we confirm that the series of the partition functions (I1.14p satisfy the Toda 
equation (11.171) . In Section 4, we provide a proof of the statement. 
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2 Quantum torus Lie algebra 

Let Vm \ where k = 0, 1,2, ■ ■ ■ and m G Z, be a set of operators which are defined by the 
following generating function for each k. 

:^{qlz)r{q~-^z): = z"™- KI'^^ (2.1) 

This formula yields the following expression of vi'K 

Vi'^ = q-^ { ^.z^ : ^{q'^zW^f-^z) : , (2.2) 
J 27r« 

where the integral means taking the residue at z = 0. This integral can be evaluated by 
plugging the mode expansion of ?/'(z), '?/'*(2;) into the generating function. Thereby, the right 
hand side of equation (12.21) is read as 

The operator (11.241) and the U{1) current (ll.20p are represented as 

= K^'^ , = Ki") . (2.4) 

Henceforth, taking the viewpoint of integrable systems, we call {k > 1) hamiltonians. 

We note that the normal ordering in the formula (12. ip is redundant when k ^ 0. Actually, 
the fermion bilinear form —ip*{q~^z)ilj{q2z) is regularized solely by the point splitting^ z — >■ 



^We take < g < 1. 
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q'^^z, without any normal ordering. Effect of the normal ordering is known from fll.2ip and 
becomes the subtraction of a finite term as 

k 

- ^Ij*{q-^z)^lj{qh) =: ilj{qh)^*{q-h) : . (2.5) 

(k) 

Thereby, the normal ordering only makes the finite gap between Vq and that obtained without 
the normal ordering as follows. 

Vo'^ = -E^'™^-^- + rz^- (2-6) 

m 

The operators Vm^ satisfy quantum torus Lie algebra (the sine-algebra [27]). The following 
commutation relations can be found out. 

"•[/(fc) T/(0 
' rn ^ ' n 

where fc, Z = 0, 1, ■ ■ ■ and m,n E Tj. The commutation relations (12.71) become the standard 
ones by shifting V^q*-'^'' — > V^*-'^'' — for k ^ 0. The hamiltonians Hk{k > 0), where Hq = 
Vq'^^ is included, generate a commutative sub-algebra of this Lie algebra. The non-negative 
modes {Jm}^^^ and the non-positive modes {J_m}^^Q of the U{1) current do as well. In 
addition to these sub-algebras, there are two more commutative sub-algebras that are generated 
respectively by {^i'^''}fc>Q and {k!;^''}^^q. All these sub-algebras are found to relate with one 
another. 

The appearance of quantum torus Lie algebra might be unexpected. However it can be 
explained from the viewpoint of the sine-algebra |27j. The sine-algebra is obtained from sl{N) 
by taking the large N limit of its trigonometrical basis. Let X, Y be two N xN unitary matrices 
given by 

X = Y,E,,+i + En,i, Y = J2^'''Ei^i, (2.8) 
1=1 1=1 

where u; is a A^-th root of the unity. These two matrices satisfy the relation XY = ujYX. Their 
non-commutative monomials X'^Y^ for < m < and < < A^ give the trigonometri- 
cal basis of sl{N). The sine-algebra is the Lie algebra of X"^Y'' at A^ = cxo and is identified 
with quantum torus Lie algebra. It is the Lie algebra derived from quantum two-torus (non- 
commutative two-torus), that is, an unital algebra with two generators U,V satisfying the 
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Im — kn l'm — kn._{(U_\l\ „ Q \ 



relation UV = qVU. Here q is regarded as the non- commutative parameter. The trigonomet- 
rical basis X"^Y'' are transmuted to the non-commutative monomials U^V^. Let us normalize 
them as follows. 

= g-¥^7'"y^ (2.9) 

These normalized ones satisfy the commutation relations (12.71) . apart from the shift of the zero 
modes. 

The sine-algebra is the algebra of trigonometric basis of s/(oo). Making use of the embedding 
sZ(oo) c gl{oo), as utilized in (12. Sp for finite A^, a trigonometric basis can be realized in terms 
of g-difference operators with respect to z. Among them, the fundamental operators are z and 
q-zd/dz _ gxp(— loggz^), which are the counterparts of X and Y . These two operators satisfy 
the relation zq~^'^^'^^ = qq~^^^'^^z. The normalized basis (12.91) of the sine-algebra correspond to 

= q-^z'^q-kz^ _ (2.10) 
The operators Vm"^ are nothing but the second-quantizations of fm"* by means of the fermions. 

Vj^^ = j^^:^{z)vt^r{z): . (2.11) 

Since we have restricted A; > 0, the operators Vm^ or the basis (I2.10p generate only the 
half of quantum torus Lie algebra. Needless to say, the operators for A; < are obtainable 
from the generating function (12.11) by choosing k as such, however those play no role in this 
article. The above half is, so to speak, a quantum cylinder Lie algebra to be obtained from a 
quantum cylinder. In the context of random plane partition, this quantum cylinder becomes a 
classical cylinder C* at the thermodynamic limit or the semi-classical limit; the Seiberg-Witten 
hyper-elliptic curves of five-dimensional A/" = 1 supersymmetric gauge theories emerge as the 
double cover of the cylinder [m [16] . See [HI HE] for the details. 



3 One-dimensional Toda chain 

In this section we show that the series of the partition functions (I1.14p provides a solution of the 
one-dimensional Toda chain. We will identify the partition functions with dynamical variables 
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bp on a Z-lattice by 



r{t;p+l] 



q 



(p+i) 



1)2 Zp+l{t) 



pez. 



r{t;p) ^ Zp{t) ' 

Then, writing x = ti, the variables (^p satisfy the following Toda equation. 



d': 



peZ. 



(3.1) 



(3.2) 



3.1 Transformations by adjoint action 

We argue transformations of the generators Vm^ by the adjoint action of the propagators G±. 
For this purpose, we conveniently start with the transformations of the fermions tlj{z),tlj*{z). 
These follow from their rotations by the transfer matrices r±(m), since G± are just the products 
of r±(m). The rotations are easily computed by recalling that each modes of the U{1) current 
rotates the fermions as follows. 



-dJrr 



(3.3) 



For instance, the rotation of ip{z) by r+(m) is computed as follows. 

+ 00 



r+(m) ipi^z) r+(m) 



+ 00 l„k„~k(rn+l) 



k=l 



k=l 

Hz) 

(1 - zg-(™+^))-V(^) • 



e * 



(3.4) 



By using this formula reiteratively according as fll.32p . we can compute the adjoint action of 
G+. In this way, we eventually obtain the following transformations of the fermions. 



G+^(z)(G+)-^= n(i-^r-^)-v(^) 



m=l 

+ 00 



(3.5) 



m=l 

+ 00 

{G.)-'Hz)G^ = n (1 - z-'q'^^h Hz) , 

m=0 

+ 00 

{G.)-^^*{z)G. = n ~ ;z-'g'"+5)-i ^*(z) . 



(3.6) 



m=0 
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Nextly we examine the transformations of the generating function (12.11) of the quantum 
torus Lie algebra. These transformations are obtained by using formulas (13.51) . (13. 6p taking 
(12. 5p into account. We illustrate the computation for the case of G+. 

k 



z{l — q'' 

k 
2 



= - n (1 - '^'^''^) r{q~'^z)i^{q^z) + . (3.7) 

z{i — q ] 

m=l ^ ' 

The last line in (13.71) can be rewritten in terms of the generating function itself. We thus obtain 

{ : i,{q^z)r{q-^^z) : "^^Y^} 

( h. 

'^-m\ u \ I*/ \ q^ 



n (1 - zq'^-n { : ^P{qh)r{q--^z) : -^^j^ \ • (3.^ 



m=l 

The similar computation goes as well for the case of G- and gives 
(G_)-^<!:^(gt.)^*(g-t.):-^^^^|G_ 



= n (1 - z-'q-'^^n : ^(gt^)^*(g-t^) : - . (3.9) 

m=l I, Z\ q ) j 

The transformations of Krf'' can be obtained from formulas ( 13.80 . ( 13.90 by reading the coef- 
ficients of the Laurent expansions of the equations around z = 0. It is evident but nevertheless 
surprising that G± generate automorphisms of the Lie algebra (12. 7p by adjoint action. Let 
us concentrate on the transformations of the hamiltonians Hj. = Vq''\ These are read from 
formulas (13. Sp . (13.90 as follows. 

G^V,^'\G^)-' = gt /^n(l-^g^-):^(gtz)^*(g-tz):, (3.10) 

m=l 
k 

(G^rX^'^G^ = JJ(l-,-V^+-):^(gl^)^*(g-|^):, (3.II) 

m=l 
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where the integrals denote taking the residue at z = 0. These integrals can be evaluated by 
using the g-binomial theorem: 

k k 



n(i+ 



xq 



E »(»+!) 
q 2 



m=l 



i=0 



k 

i 



X 



(3.12) 



where 



k 

i 



q {q : q)i{q : q)k- 



(g:g)„ = (l-g)(l-g2)...(i_g"). 



(3.13) 



By expanding the products in the right hand sides of equations f l3.10l) . fl3.11l) . these are brought 
to the following form. 



i(k — i) 



i=0 
k 

E(-)V 

i=0 



k 

i 



(k) 



i{k — i) 



k 



V 



(k) 



(3.14) 
(3.15) 



3.2 Toda equation 

The transformations of Hi = Vq^'^ in formulas (13. 14p . ( 13.151) are read as 



(3.16) 
(3.17) 



These transformations deal naturally in the evolution of the partition function (11. 9p by the time 
X = ti. The operator G^e^^^^G^ that appears in the expression (ll.3ip evolves according to 



It can be also written as 



These two descriptions lead to 



{p\{V^'^-V}'^)G^e^('^G^ \p) 



(3.18) 



(3.19) 



(3.20) 
(3.21) 
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Let us derive the Toda equation fl3.2p . Owing to the identification (13. II) . it suffices to prove 
the following identity. 

dlZj, - [d^Z^f = Vi Zp-i, pel. (3.22) 

We first rewrite the left hand side of equation ( 13.22^ , using the expression ( 11.3ip and applying 
formulas fl3.18l) . fl3.19p in it, as follows. 

^p^x^p ~ (dxZpY 

= Z, . {p\ iV^'^ - Vl'^) G+e^WG_ iV^'^ - V^\^) \p) 
- {P\ iV,^'^ - V}'^) G+e^WG_ \p) ■ {p\ G+e^WG_ {V,^'^ - y«) \p) (3.23) 

The matrix elements in the right hand side of this equation can be translated to the fermion 
correlation functions, by replacing V^''^'' — V^i with the states generated by these operators. 
The corresponding states are 

(Vji) _ v«) ^ + q^^-h-p-ir^ \P) (3.24) 

and its conjugate state. Thereby, we can eventually translate the right hand side of equation 
(13.231) into the following combination of the correlation functions. 

^P ^x^P ~ (dxZp)'^ 

= q'^+'Zp ■ {p\ G+e^(*)G_ \p) 

- g^P+^pl ^.pr^+i G+e^WG_ \p) ■ {p\ G+e^^G^ \p) (3.25) 

Wick's theorem shows that correlation functions of free fermions are factorized into products 
of their two point functions. The four point function in the right hand side of equation fl3.25p 
is factorized into 

^{p\ i^-p%+i G+e^WG_ \p) 

= ^{P\ G+e^WG_ \p) ■ ^{p\ G+e^WG_ \p) 

Zjp Zjp 

+^{p\ G+e^WG_ \p) ■ ^{p\ G+e^WG_ \p) . (3.26) 

Zjp Zjp 
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In the right hand side of this equation, the first term equals, making use of the relations 
|p + 1) = ^/'_p_i|p) and \p — 1) = ipplp), to Zp^Zp^iZp^i. Thus, Wick's theorem leads to 

Zp ■ {p\ i^-pr,+i G+e^(*)G_ \p) 

= Zp+^Zp^, + {p\ i^-p%+i G+e^WG_ \p) ■ {p\ G+e^(*)G_ V'-p-i^; \p) • (3.27) 

By plugging this formula into the right hand side of equation (13.251) . we obtain q^'^^^ Zp+iZp^i. 
Thereby this completes the proof. 



4 One-dimensional Toda hierarchy 

In this section we prove that the series of the partition functions fll.l4p is a tau function of the 
one-dimensional Toda hierarchy. 

We first recall the theory of tau functions of the Toda hierarchy [28l [29| [30] . The two- 
dimensional Toda hierarchy has two series of commuting flows and thereby two series of time 
variables, T = (Ti, T2, ■ ■ ■ ) and T = (Ti, T2, ■ ■ ■ ), each of which describes each the commut- 
ing flows. Tau functions of the two-dimensional Toda hierarchy are admitted to have several 
expressions including the realization by means of free fermions or free bosons. The standard 
description in terms of free fermions is 

where g is an element of GL{oo). Ck and are numerical constants which originate in the ambi- 
guity of the tau function. The two-dimensional Toda hierarchy reduces to the one-dimensional 
Toda hierarchy when the two-sided time evolutions degenerate. This reduction imposes the 
following constraint on the tau function. 

(A + A)^""*^-^^"'""- " = 

When the condition fulfilled, the tau function translates to the tau function of the one- 
dimensional Toda hierarchy. Owing to the degeneration, the one-dimensional Toda hierarchy 
has one series of commuting flows. The corresponding time variables can be identified with 
T = (Ti, T2, ■ ■ ■ )■ The tau function has the following expression in terms of free fermions. 
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The coupling constants t = {ti,t2, ■ ■ ■) in the series of the partition functions (11.141) is eventually 
identified with the standard Toda time variables T = (Ti, T2, ■ ■ ■ ) by = {—Ytk. 
To implement the constraint fl4.2p . g is chosen to satisfy the constraint 

Jkg = gJ-k, k = l,2,---. (4.4) 
Under this condition, the foregoing expression of T^'^°'^°-{T\p) can be rewritten as 

This shows that ri™"(T;p) is a tau function of the modified KP hierarchy as well. 
4.1 Shift symmetry 

Among automorphisms of the Lie algebra (12. 7p . we pay special attention to the following shift 
symmetry: 

k k 

and Vm^ are left unchanged. Using this symmetry, three commutative sub-algebras generated 

respectively by {K)^^''}fc>o' {^fc'^^}fc>o {^-fc^}fc>o become conjugate to one another. 

The symmetry (14.61) becomes eventually one of the automorphisms generated by the adjoint 
action of G±. Combining the transformations (13. 81) . (13. 91) . we find 

= H'z' |: : "7(13^} • (4-7) 

Taking account of the mode expansion (12.11) . we can read the symmetry (14. 6 p from formula 
(14. 7p in the following form. 

This formula shows that the conjugacy between the three sub-algebras is realized as their 
transformations by (G'_G'+)^^. We shall describe such transformations in some detail. Let 
k > 1. Putting m = in (14.80 . we obtain 

(G_G+) (v,^'^ - ^) (G_G+)-i = {-)'V^'^ . (4.9) 



1 - 
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Similarly, putting m = —k in (14.81) . we obtain 

(G-G^)-' (v^'^ - = (-) Vit^ . (4.10) 

Let us consider two more commutative sub-algebras which are generated respectively by 
{^i°^}fc>o ^^"^ {^-A:^}fe>o' These sub-algebras become conjugate to the aforementioned three 
sub-algebras by the adjoint actions of (G'_G'+)'^^ and . To see this, note that W (I1.35P is 
able to rotate vjj} to V^j^ by the adjoint action. Actually, we have the following transforma- 
tions. 

,fvf),-f = If, (4.11) 
,-fV'W,¥ = (4.12) 

W 

To obtain these formulas, note that, by the adjoint action of g^, the fermions transform as 
q 2 Wraq 2 = g 2 -^p^ and q 2 2 = g 2 7/;^. By usmg these transformations, the left hand 

side of (14. lip can be computed as 

q^ q ^ = g 2 2^ g 2 : ?/^fc_„V^„ : q 2 
= g ^ }_^g'="+^ -2 : ^,_„z^; : 

neZ 

= (4.13) 

which is nothing but V^^\ Thus we obtain the formula (14. lip . The similar computation leads 
to the formula (I4.12p . 

The formulas (14. lip . (I4.12p in addition to fj49|) . (I4.10p show that all the five sub-algebras 
are conjugate to one another by the adjoint actions of (G'_G'+)^^ and . Therefore, the 
hamiltonian i/^. = V^*''^'' can transform to J±k = V^/j by the adjoint actions of (G'_G'+)^^ 
and g^^. Actually, as seen from formulas (14.90 . (14. lip . g^(G_G'_|_) rotates to J^, while 
q~^ (G^G^)^^ rotates iffc to J_fc, as seen from formulas (I4.10p . (I4.12p . We write down these 
transformations in the following form convenient for the later use. 

G^ (vi'' - (G^r = {-nq'^G^yKj.iq^G^), (4.14) 

(G-)-^ (^^'^ - Y^) = (-)^'(G+g^)J_.(G+g^)-^ (4.15) 
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4.2 The proof 

4.2.1 New representation of the partition functions 

We derive the expression ( ll.lSp of the partition function ( 11.91) . 
Let us first rewrite G^^e^^^^G^ as follows. 

= G'+e5^(*)(G'+)-^G+G_ (G.j-^e^^WG.. (4.16) 

The transformations of e^^*-*-* by the adjoint actions of G+,(G_)~^ in this expression can be 
evaluated by using formulas fl4.14p . (14.151) . Eventually, these transformations are expressed as 

G+e^^W(G+)-i = e^^"^-^ (gfG_)-ie^Er^i(-)'*^^^ (gfG_) (4.17) 
(G_)-ie5^WG_ = e^^-^ (G+gf)e^S-^(-)'^*^^- (G+g^?)-^ (4.18) 
By plugging these expressions into the right hand side of (14.161) . we obtain the following formula. 

G+e^«G_ = e^^-^ (g^G_)-i e5J^-i(-)'=*^^^ 

xg^ e^^^-(-)'*^^- (G+g^)-\ (4.19) 

where is the element of GL(oo) given by (11.34p . 

Making use of this formula we arrange the expression (11.310 as 

Z,{t) = (p|G+e^WG_b) 

^ gEr^i^ (G_)-ig-f e^S^"-^-)'*^^^ g, ei^^-i(-)'*'^-^- g'f (G+)-i \p) 

y^oo _tfe9 1 ^oo /- ^fcJ- 7 1 v^oo I \k-i J W , , / , X 

= e'^'==^^ e^Efe^i^-'' g^ e2^'==i(~) q~~ \p) . (4.20) 

Considering in the last line that the W^-charge of the state \p) is \p{p + l)(2p +1), we finally 
obtain the expression (ll.lSp . 

4.2.2 Reduction to one-dimensional Toda hierarchy 

Let us show that T(t;p) = q^^^P^^^^'^^'^^^ Zp(t) is a tau function of the one-dimensional Toda 
hierarchy. The key is the fact that g^, satisfies (14. 4p . This can be seen by using formulas (14. 9p . 
flCTD . fimi) . fHl2D as follows. 

Jfcg* = Jkq~iG-G+yq- 



w 

2 
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(k) 




W W 

q-{G^G+fq-J^k 



gi. J-k ■ 



(4.21) 



The constraint (14.41) is equivalent to the constraint (14.21) on the tau function given by (14. ip 
taking = g^. This means that actually gives a solution of the one- dimensional Toda 
hierarchy. Therefore it follows from (14.31) that T{t;p) is the corresponding tau function by the 
identifications = {—Ytk- This completes the proof. 

5 Conclusion and discussion 

We investigated melting crystal, which is known as random plane partition, from the viewpoint 
of integrable systems. We proved that a series of partition functions of melting crystals gives 
rise to a tau function of the one-dimensional Toda hierarchy, where the models are defined by 
adding suitable potentials, endowed with a series of coupling constants, to the standard statis- 
tical weight. We showed that these potentials are converted to a commutative sub-algebra of 
quantum torus Lie algebra. Further exploiting the underlying algebraic structure, a remark- 
able connection between random plane partition and quantum torus Lie algebra was revealed. 
This connection substantially enabled to prove the statement. Based on the result, we briefly 
argued the integrable structures of five-dimensional M = 1 supersymmetric gauge theories and 
A-model topological strings. The aforementioned potentials correspond to gauge theory ob- 
servables analogous to the Wilson loops, and thereby the partition functions are translated in 
the gauge theory to generating functions of their correlators. In topological strings, we particu- 
larly comment on a possibility of topology change caused by condensation of these observables, 
giving a simple example. 

In four-dimensional M = 2 supersymmetric gauge theories, the authors of [10] obtained the 
generating function of correlation functions of the higher Casimir operators in the fermionic 
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form 

= (p|es-^ie^*°°=o (^^'=+^6^-^-1 (5.1) 

where Vk are the fermion bihnear forms introduced by Okounkov and Pandharipande [23], 
and Xk are the couphng constants of the higher Casimirs of the gauge theory. The above 
partition function also appears in the Gromov-Witten theory as the generating function of the 
absolute Gromov-Witten invariants on CP^ [23l[T8]. Before this fermionic representation was 
presented, Getzler had conjectured, and later proven [21], that the generating function is a tau 
function of the one-dimensional Toda hierarchy. Getzler's proof is, however, fairly complicated 
and somewhat indirect, combining the Virasoro conjecture \2A\ with the partial result that the 
Toda equation holds on the subspace X2 = = ■ ■ ■ = 0. It will be therefore an interesting 
problem to give a more direct proof on the basis of the fermionic representation. A possible 
scenario will be, as we have done in the five- dimensional case, to find a suitable analogue of 
and to rewrite the foregoing fermionic representation into a standard form like (14.31) and (14. 5p . 
Unfortunately, as we remarked in the end of section 1.5.2, the naive four- dimensional {R 0) 
limit of itself does not work. Since the role of g* reminds us of various "dressing operators" 
in the work of Okounkov and Pandharipande [231 [25] , a correct four-dimensional analogue of 
g^ might be hidden therein. 
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